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Abstract. LISA is a joint space mission of the NASA and the ESA for detecting low 
frequency gravitational waves in the band 10 -5 — 1 Hz. In order to attain the requisite 
sensitivity for LISA, the laser frequency noise must be suppressed below the other secondary 
noises such as the optical path noise, acceleration noise etc. This is achieved by combining 
time-delayed data for which precise knowledge of time-delays is required. The gravitational 
field, mainly that of the Sun and the motion of LISA affect the time-delays and the optical 
links. Further, the effect of the gravitational field of the Earth on the orbits of spacecraft is 
included. This leads to additional flexing over and above that of the Sun. We have written a 
numerical code which computes the optical links, that is, the time-delays with great accuracy 
~ 10~ 2 metres - more than what is required for time delay interferometry (TDI) - for most 
of the orbit and with sufficient accuracy within ~ 10 metres for an integrated time window of 
about six days, when one of the arms tends to be tangent to the orbit. Our analysis of the 
optical links is fully general relativistic and the numerical code takes into account effects such 
as the Sagnac, Shapiro delay, etc.. We show that with the deemed parameters in the design 
of LISA, there are symmetries inherent in the configuration of LISA and in the physics, which 
may be used effectively to suppress the residual laser noise in the modified first generation 
TDI. We demonstrate our results for some important TDI variables. 
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1. Introduction 

A number of ground-based large-scale interferometric gravitational wave detectors, with optimal 
sensitivity in the frequency window ~ 10 Hz - 1 kHz are operational world-wide pp. Among the 
large scale detectors, detectors of armlengths of few kilometres, the LIGO detectors of the US 
have now been continuously operating at initial design sensitivity which gives them a maximum 
range of about 25 Mpc for compact neutron star - neutron star binaries. The VIRGO detector 
of France and Italy has also attained comparable sensitivity. The detectors have been built with 
the realistic goal of directly observing gravitational waves (GW's) for the first time. 

A natural limit occurs on decreasing the lower frequency cut-off of 10 Hz because it is not 
practical to increase the armlengths on ground and also because of the gravity gradient noise 
which is difficult to eliminate below 10 Hz. The solution is to build an interferometer in space, 
where such noises will be absent and allow the detection of GW in the low frequency regime. 
LISA - Laser Interferometric Space Antenna - is a proposed mission which will use coherent laser 
beams exchanged between three identical spacecraft forming a giant (almost) equilateral triangle 
of side 5 x 10 6 kilometres to observe and detect low frequency cosmic GW [2]. The ground- 
based detectors and LISA complement each other in the observation of GW in an essential way, 
analogous to the optical, radio, X-ray, 7-ray etc., observations do for the electromagnetic waves. 

In ground based detectors the arms are as symmetrical as possible so that the laser light 
experiences nearly identical delay in each arm of the interferometer. This arrangement reduces 
the laser frequency /phase noise at the photodetector. Reduction of noise is crucial since the raw 
laser noise is orders of magnitude larger than other noises in the interferometer. But perfect 
symmetry is not possible, and an efficient system of servo loops is necessary for reaching a 
noise level compatible with the required sensitivity. The required sensitivity of the instrument 
can thus only be achieved by near exact cancellation of the laser frequency noise plus a good 
symmetry of the arms. However, in LISA, the lack of symmetry will be much larger than in 
terrestrial instruments, and the laser noise, though reduced by stabilisation techniques (still to be 
demonstrated), will probably be still too high. LISA consists of three correlated interferometers, 
which produce redundancy in the data, and this can be used to suppress the laser frequency 
noise. In LISA, six data streams arise from the exchange of laser beams between the three 
spacecraft - it is not possible to bounce laser beams between different spacecraft, as is done in 
ground based detectors, because after 5 million km propagation, the intensity of light reaching 
the target spacecraft is reduced by 10 orders of magnitude; but in the target spacecraft a laser 
is locked in phase on the received wave, so that the secondary beam is re-emitted without loss 
of phase information to the primary source. This is analogous to the RF transponder scheme, 
as was done in the early experiments for detecting GW by Doppler tracking a spacecraft from 
Earth [3]. 

Laser frequency noise which dominates the other noises by 7 or 8 orders of magnitude 
must be removed if LISA is to achieve the required sensitivity of h ~ 10~ 22 , where h is the 
metric perturbation caused by a gravitational wave. This cancellation is achieved by time-delay 
interferometry (TDI) where the six data streams are combined with appropriate time-delays. 
This is possible because of the redundancy present in the data. This work was put on a sound 
footing by showing the data combinations had an algebraic structure; the data combinations 
cancelling laser frequency noise formed the module of syzygies over the polynomial ring of time- 
delay operators [I] . This work was done for stationary LISA in flat spacetime where the motion of 
LISA as well as the ambient gravitation field, mainly that of the sun, was ignored. These were the 
so-called first generation TDI. However, LISA spacecraft execute a rotational motion and also the 
background spacetime is curved, all of which affect the optical links and the time-delays. Thus 
the Sagnac effect, Einstein effect, Shapiro delay, etc. are important and must be incorporated 
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into the analysis if the laser frequency noise is to be effectively cancelled. We here take into 
consideration all these effects in the full framework of general relativity. However, we compute 
the orbits of spacecraft in the Newtonian framework. The base orbits we take to be Keplerian in 
the gravitational field of the Sun only, assuming the Sun to be a point mass. On these base orbits, 
we linearly superpose the perturbative effect of the Earth's gravitational field. We choose the 
Earth over Jupiter firstly because, the Earth perturbs the Keplerian orbit in resonance, resulting 
in unbounded growing of the perturbations and secondly, because of the technical reason that 
the Earth's effect can be easily accommodated within Clohessy- Wiltshire (CW) framework [5]. 
Moreover we argue that Jupiter's effect is less than 10% of that of the Earth's and hence not a 
dominant one. The analytic approach helps to gain insight and understanding of the problem. 

Finally, given the arm flexing for our model of LISA, we compute the residual laser frequency 
noise spectrum for some important TDI observables, namely, the Sagnac, the Michelson and 
the Symmetric Sagnac in their modified first generation form. We find that the residual laser 
frequency noise in general tends not to be very high as compared with the secondary noises. If 
this level of noise is found to be acceptable, then there may be no need to use second generation 
TDI observables, which in general involve higher degree polynomials in time-delay operators and 
thus require more interpolations which in turn result in larger errors in the data analysis. 

We believe that the computations that we present here of the model would be of help in the 
development of a LISA simulator, the LISACode for instance [6], because, (i) we have taken into 
account the relativisitic effects and (ii) the effect of the perturbation of the orbit due to the Earth. 
Also it would be useful to compare the model with actual data and look for any discrepancies. 
Any discrepancy arising could be interesting because that would imply the existence of some 
physical cause which has been overlooked and therefore would have to be incorporated into the 
data analysis. 



2. The spacecraft orbits and flexing of LISA's arms 

We first describe the orbits in the gravitational field of the sun only. These are the usual Keplerian 
orbits. We then give the description of the same orbits in terms of the CW equations. This paves 
the way for including the effect of the Earth. Finally, we use the CW framework to include the 
perturbative effects of the Earth. 



2.1. The Keplerian orbits of spacecraft in the Sun's field 

The Keplerian orbits are chosen so that the peak to peak variation in armlengths is the least 
~ 48000 km, see [?]■ We summarise the results of this paper below. We choose the Sun as the 
origin with Cartesian coordinates {X, Y, Z} as follows: The ecliptic plane is the X — Y plane 
and we consider a circular reference orbit of radius R equal to 1 A. U. centred at the Sun. Let 
5q = 5a/8 where a = Lq/2R and Lq ~ 5,000,000 km is a constant representing the nominal 
distance between two spacecraft of the LISA configuration. We choose the tilt of the plane of 
the LISA triangle to be S = ir/3 + S (this results in minimum flexing of the arms). We choose 
spacecraft 1 to be at its lowest point (maximum negative Z) at t = 0. This means that at this 
point, Y = and X ~ R(l — e). The orbit of the first spacecraft is an ellipse with inclination 
angle eo, eccentricity e and satisfying the above initial condition. 
From the geometry, eo and e are obtained as functions of S, 

a sin S 
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The equations for the orbit of spacecraft 1 are given by: 

X\ = i?(cosi/>i — e)coseo, 
Yi = Ryl — e 2 sin-0i, 

Z\ = — i?(cos^i — e) sine - (2) 
The eccentric anomaly ip\ is implicitly given in terms of t by, 

■01 — e sin ipi = fit — (j)Q , (3) 

where t is the time and ft is the average angular velocity and <f>o the initial phase. The orbits 
of the spacecraft 2 and 3 are obtained by rotating the orbit of spacecraft 1 by 27r/3 and 4-7r/3 
about the Z— axis; the phases ^^z-, however, must be adjusted so that the spacecraft are at a 
distance ~ Lq from each other. The orbital equations of spacecraft k = 2, 3 are: 

X k = Xi cos a k - Y\ sin a k , 
Yfc = X\ sin cTfc + Yi cos a k , 

Z k =Z x , (4) 

where a k = (k — 1) with the caveat that the ipi is replaced by the phases ip k , where they are 
implicitly given by, 

ip k - e sin tp k = fit - cr fc - <j> . (5) 

These are the exact (Keplerian) expressions for the orbits of the three spacecraft in the Sun's 
field. In [7] it was shown that these orbits produce minimum flexing of LISA's arms when only 
the Sun's field is considered. Our next goal is to include the Earth's field and compute the 
flexing of LISA's arms in the combined field of the Sun and Earth. For this purpose we use the 
Clohessy-Wiltshire framework. 



2.2. The Clohessy-Wiltshire framework 

Clohessy and Wiltshire make a transformation to a frame - the CW frame {x, y, z} which has its 
origin on the reference orbit and also rotates with angular velocity fl. The x direction is normal 
and coplanar with the reference orbit, the y direction is tangential and comoving, and the z 
direction is chosen orthogonal to the orbital plane. They write down the linearised dynamical 
equations for test-particles in the neighbourhood of a reference particle (such as the Earth). 
The length scale here is the Earth-Sun distance of 1 A. U. and the equations are applicable to 
distances small compared with this length scale. Since the frame is noninertial, Coriolis and 
centrifugal forces appear in addition to the tidal forces. The advantage of the CW equations is 
that it is easy to see that to the first order in a (or equivalently e) there exist configurations 
of spacecraft so that the mutual distances between them remain constant in time. The flexing 
appears only when we consider second and higher order terms in a. In fact in [7] we find that 
the second order terms are sufficient to describe the flexing of LISA's arms quite accurately. 

We take the reference particle to be orbiting in a circle of radius R with constant angular 
velocity Q. Then the transformation to the CW frame {x,y,z} from the barycentric frame 
{X, Y, Z} is given by, 

x = (X - R cos fit) cos fit + (Y - R sin fit) sin fit , 
y = - (X - R cos fit) sin fit + (Y — R sin fit) cos fit , 

z = Z. (6) 
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In the CW frame, the CW equations to second order in a (this includes upto the octupole field 
of the Sun) are, 

x - 2Qy - 3Ci 2 x + ^-(2x 2 - y 2 - z 2 ) = , 
■■ 6aft 2 

y + 2ilx — xy = , 

Lq 

z + tt 2 z — xz =0. (7) 

Lq 

If we drop the terms in a in these equations we get the original C W equations (upto quadrupole) . 
The solutions to these equations (that is upto quadrupole order) we call the zero'th order. Among 
these we choose the solutions which form an equilateral triangular configuration of side Lq. For 
the fcth spacecraft we have the following coordinates: 

x k = - ip cos(fW -a k - M , 
y k = po sin(Ot - a k - 4> ) , 

Zk = — y Pq cos (^ f -<*k- 4>o) , (8) 

where po = Lo/y/3. Also at t = the initial phase of the configuration is described through 4>q. 
In this solution, any pair of spacecraft maintain the constant distance Lq between each other. 

In [7] we have shown that if include the a terms (octupolar terms) and solve perturbatively 
using the zeroth order solution as given by Eq.®, we obtain the flexing of the arms. Further, 
this approximate solution agrees to a remarkable degree with the flexing deduced from the exact 
Keplerian orbits. 



2.3. The effect of the Earth 

LISA follows the Earth 20° behind. We consider the model where the centre of the Earth leads the 
origin of the CW frame by 20° - thus in our model, the 'Earth' or the centre of force representing 
the Earth, follows the circular reference orbit of radius 1 A. U. Also the Earth is at a fixed position 
vector r® = (a;©,y®,z©) in the CW frame. We find that a;© = —R(l — cos 20°) ~-9x 10 6 km, 
y® = i?sin20° ~ 5.13 x 10 7 km and z® = 0, where we have taken R to be 1 A. U.. The force 
field F due to the Earth at any point r (in particular at any spacecraft) in the CW frame is given 
by: 

F(r) = -GMe^^3, ( 9 ) 

where M® ~ 5.97 x 10 24 kg is the mass of the Earth and G = 6.67 x 10~ n kg _1 m 3 sec~ 2 Newton's 
gravitational constant. 

In order to write the CW equations in a convenient form we first define the small parameter 
e in terms of the quantity w® = GM®/d®, where d® = |r®| is the distance of the Earth from the 
origin of the CW frame; d® ~ 5.2 x 10 7 km which is more than 50 million km. So when deriving 
the forcing term we make the aprroximation |r — r®| w rf®, that is, we neglect |r| compared to 
rf®. It will turn out that the flexing due to the Earth is small so that this approximation is not 
unjustified. We define e = u^/fl 2 ~ 7.16 x 1CP 5 which is the just the ratio of the tidal forces 
due to the Earth and the Sun. The CW equations including the Earth's field take the form: 

x - 2Qy - 3Q 2 x + eVl 2 (x - x m ) = , 

y + 2nx + en 2 (y- y(B ) =0, 

z + n 2 {l + e)z =0. (10) 
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Note that the compounded flexing due to the combined field of Earth and Sun is nonlinear; it 
is infact a three body problem. We however solve this problem approximately. Assuming that 
both effects are small we may linearly add the flexing vectors due to the Sun and Earth; that 
is, add the perturbative solutions obtained from Eqs.10 and ifTOlk the nonlinearities appear at 
higher orders in a and e. These would modify the flexing but we may neglect this effect because 
of the smallness. As it will turn out, the flexing produced due to the Earth is of the order of 1 
or 2 m/sec upto the third year, just about 40 % of that due to the Sun. But, as shown in [7] the 
flexing produced by the Sun's octupole field is nearly exact to that produced by the Keplerian 
orbits. Thus we may do better by linearly adding the flexing vector produced by the Earth to the 
Keplerian orbit of the relevant spacecraft. We therefore, first compute the motion of spacecraft 
perturbatively using the zero'th order solutions as given in Eq.®. This will induce a flexing of 
the LISA arms only by the Earth's field. 

We now seek perturbative solutions to Eq. lfT0|) to the first order in e. We write, 
x = xq + ex\, y = yo + eyi, z — zq + ezi where Xo,yo, zq are solutions at the zeroth order given 
by Eq.|(8]). We put = (or equivalently include it in 0o) in these solutions for simplifying the 
algebra. 

Note that the z equation is decoupled from the x and y equations which are themselves 
coupled; infact the z equation can be solved exactly. We also assume the initial conditions for 
the perturbative solutions to be homogeneous, that is, we take, x\ — y% = Z\ = x\ = iji = z\ = 
at t = - the spacecraft are in the desired positions initially. 

We first solve the x and y equations. To the first order in e, the equations for these 
perturbations are: 

xi — 2ilyi — 3tt 2 xi = fl 2 XQj + — fl 2 po cos(ili — 0o) , 

y\ + 2Qxi = Sl 2 y© - f2 2 po sin(f2i - O ) . (11) 

We note that the forcing terms on the right hand side of these equations appear at the same 

frequency SI and hence they imply resonance. This means that the Earth's effect on LISA is 

cumulative as the detailed calculations show below. Therefore, it is most important to include 

the effect of the Earth on LISA. 

The equation for y\ can be easily integrated, once, with the initial conditions mentioned 

above and the y\ substituted in the x\ equation resulting in a decoupled equation for x\. This 

decoupled equation can be solved to yield the solution for the perturbation x\ . The solution x\ 

in turn can be substituted back into the y\ equation to obtain a first order equation for y\ and 

integrated with the initial conditions. Without further ado we state the results: 

5 

x\ = — Pe cos(Qt — + X0 + 2z/0$l£ — 2pocos</>o + -po^lt sin(f2£ — </>o) , 

3 

t/i = 2 ( o e [sin(Sl£ - <j> @ ) + sin0 e ] - -p [sin(f2t - <j> ) + sin</> ] 

5 3 
+ -p Qt cos(Qt - (j) ) - fM(2iE ffi - 3p cos fa) - - fl 2 t 2 y e , (12) 

where, 

g 

P% = (x@ - 2p cos O ) 2 + (2y e - ~^Po sin (j) ) 2 , 
2yg - fpo sin 0o 

tan0 e = ^ — . 13 

- 2p cos 0o 

The z equation can be exactly integrated and used directly to obtain the flexing. However, we 
can also expand this solution to the first order in e and the result is: 

z\ = po [ilt sin fit cos 0o — (fltcosilt — sinftt) sin0o] . (14) 
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We observe that since is very large compared to other distances in p®, we have /% ~ 2y® ~ 10 
km. Secondly 4>§ ~ ir/2 for the same reason. For spacecraft 1, for instance, with the initial 
condition O = 0, 0® ~ 91.83°. This aids in simplifying much of the computations. 

We now turn to the flexing of the arms for which we must compute the perturbation of each 
spacecraft orbit due to the Earth. As argued previously, we merely add the perturbation given 
by eri = e(xi,yi, z%) to the Keplerian orbit of each spacecraft. For spacecraft 1, the Keplerian 
orbit is given by Eq. ((5J) and Eq. ([3]) in the barycentric frame. We denote this zeroth order orbit 
by the vector trajectory Rn(i; 4>o) having the initial phase 4>q, We then choose the zeroth order 
solution in the CW frame corresponding to this orbit which is given by Eq. with k = 1 or 
<7i = 0. With this solution we compute the perturbative solution ri(i; 4>o) and finally obtain the 
total orbit Ri(f;0o) = Ro(*j0o) + er i(i;0o)- We repeat the same procedure for the other two 
spacecraft and obtain the orbits which now include the effect of the Earth as well. The flexing 
of the LISA arms is now due to both the Sun and the Earth. Below in Figure [TJ we plot the 
flexing of LISA's arms as a function of t, assuming constant length L in the Sun's field, that 
is, we take Ro as given from Eq.® so in effect we are considering only the flexing due to the 
Earth. We also take 0o = as the initial condition. 
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Figure 1. The figure shows the effect of only the Earth on flexing of the three arms of LISA 
for a period of three years. The flexing grows to a maximum of about 20,000 km in the third 
year which is about 40% of the flexing due to the Sun. 

We notice that the flexing increases with time; in third year it grows to a maximum of about 
20,000 km from the initial value. This variation is about 40% of that due to the Sun - ~ 48, 000 
km. Of course, one must consider the vector additions, so that the total flexing will in general 
be less than this or even reduce if the flexing vectors are oppositely oriented. Also by looking at 
the slope of the curves we note that the rate of change of armlength is of the order of a metre/sec 
in the the second year and less than 2 metres/sec in the third year. Thus at least in the first 
three years we do not expect the flexing to be affected too much by the Earth. The flexing also 
depends on the initial phase 4>q, that is, the epoch at which the gravitational field of the Earth 
is 'switched on'. Here there is also a symmetry; the flexing profiles remain invariant if we change 
0o to 0o + 7r/3. Although, for different </>o the flexing profiles differ in detail, they display the 
similar qualitative behaviour. 

The next task is to compute the optical links with the above spacecraft orbits from which 
the total flexing due to the Sun and Earth can be deduced. 
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3. The numerical computation of optical links 

The time-delay that is required for the TDI operators needs to be known very accurately - 
at least to 1 part in 10 s , that is, to about few metres - for the laser frequency noise to be 
suppressed. In order to guarantee such level of accuracy, we have to numerically compute the 
optical links or the time-delay. This approach is guaranteed to give the desired accuracy or 
even better accuracy than what is required. In this approach, we numerically integrate the null 
geodesies followed by the laser ray emitted by one spacecraft and received by the other. This 
computation is performed in the barycentric frame, and taking into account the fact that the 
spacetime is curved by the Sun's mass. The computation here is further complicated by the 
fact that the spacecraft are moving in this frame of reference and the photon emitted from one 
spacecraft must be received by the other spacecraft. We use the Runga-Kutta numerical scheme 
to integrate the differential equations describing the null geodesies. But since the end point of 
the photon trajectory is not known apriori, an iterative scheme must be devised for adjusting 
the parameters of the null geodesic, in order that the worldlines of the photon and the receiving 
spacecraft intersect. We have devised such a scheme based on the difference vector between the 
photon position vector and receiving spacecraft position vector. The six optical links have 
thus been numerically computed with sufficient accuracy required for TDI: for 98% of the time 
the code we have devised gives excellent results to the accuracy of 10 -2 metres with 10 5 steps. 
For the rest of the time, 2%, when the any one of the arms tends to lie tangent to the orbit, one 
must increase the number of steps. We increase the number of steps to 10 7 . Then the code gives 
results accurate upto 10 metres except in a window of about half an hour when the error exceeds 
this value and becomes unacceptably large. This is because the differential equations describing 
the null geodesies encounter sign changes in the components of the tangent vector which must 
be carefully incorporated into the integration scheme. Such windows occur six times in a year 
two months apart. In this paper we choose to ignore these time windows. 



3.1. Optical links: integrating the null geodesic equations 

We now turn to the optical links. The spacetime geometry taking only the Sun as the gravitating 
mass is given by the Schwarzschild spacetime whose metric in isotropic coordinates is described 
by: 

ds 2 = f(r)c 2 dt 2 - g{r)[dr 2 + r 2 (d8 2 + sin 2 6d^ 2 )} , (15) 
where the functions f(r) and g(r) are given by: 

/(r) = l-^, g(r) = l + ^. (16) 

Here m = GM/c 2 where M is the mass of the Sun, c the speed of light and G the Newton's 
gravitational constant. 

The null geodesies satisfy the differential equations: 

(17) 

(18) 
(19) 

g r 2 sin 2 C 

Here e rj g are the appropriate signs for the orbit equations, the overdot is d/cdt, where c is the 
speed of light, b is the impact parameter and L the azimuthal angular momentum. 
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Since the spacecraft are moving, the parameters b and L are not known apriori, but should 
be obtained from the solution by an iterative scheme. Let us first consider the optical link, 
say link 12, from S/C 1 (position vector r±(t)) to S/C 2 (position vector r 2 (t)) where t is the 
coordinate time in the isotropic coordinates. Since the space is almost flat, we take the zero-th 
order estimate as flat spacetime. Then the initial estimate of b is given by: 

6(0) = |ri(t ) x n| (20) 
where to is the time at which S/C 1 emits the photon. Similarly for L, the zero-th order estimate 
denoted by L/m is the z-component of the vector r±(to) x n( ) and n( ) points towards S/C 2. 
To start with a better estimate of nm) we take the better estimate of the delayed position of 
S/C 2 at t + Lq/c (note this is still not the correct final position of S/C 2, when the photon 
meets S/C 2). The null geodesic (photon orbit) is then integrated till the time as required in 
flat spacetime. The photon obviously does not hit the S/C 2, because (i) S/C 2 has moved in 
the meanwhile, (ii) Shapiro delay: the photon is delayed because of the gravitational field of 
the Sun. The difference vector between the photon and S/C 2 drives the iterative scheme: we 
decompose this vector into parallel and perpendicular components with respect to ri( ) and use 
these projections to obtain a modified direction say nm (perpendicular component) and a new 
time of flight (parallel component). The nm produces new values of b and L, say, 6m, Lm and 
so on until convergence is reached to the desired accuracy. 

We divide our discussion into two parts: (i) normal time epochs and (ii) anomalous time 
epochs (b ~ L). The anomalous time epochs, occur sixth of an year apart when any one of the 
arms lies tangential to the orbit, that is, for the relevant null geodesic b ~ L. With the initial 
conditions, we have chosen, the first such epoch for the link 12 occurs at ~ 52 x 10 5 sees. In the 
normal case, we use 10 5 steps in the Runga-Kutta scheme, each step of about 50 km. Just 2 or 
3 iterations suffice to produce the necessary convergence and with an accuracy to 10 -2 metres. 
In the anomalous case, we need more steps in the Runga-Kutta scheme; we use 10 7 steps which 
gives an accuracy within 10 metres except for a window of about half hour. In the anomalous 
case at least one of the 8 and r (or equivalently e r and eg) change sign along the null geodesic 
and this sign change must be taken into account in the Runga-Kutta integration scheme as close 
as possible to the turning point. We achieve this by reducing the step size to ~ 0.5 km, so that 
the error remains at an acceptable level. Even then in a window of about half an hour, every six 
months, the error becomes large. We ignore these points by smoothly interpolating. 

The Figure [2] shows the comparison of the optical link L 12 just due to the Sun and when 
the Earth's field is also included. The flexing increases by ~ 20, 000 km (as compared to L Q ) 
when the Earth's field is included. The figure is drawn for a time period of about three years and 
0o is chosen to be zero. The important difference due to the combined field is that the motion 
of spacecraft and hence the variations in armlengths is no more periodic and infact grows with 
time. For certain links (in our case links 23 and 32), the flexing is infact reduced because the 
flexing vectors due to the Sun and Earth are oppositely oriented. 

Figure [3] shows all the six optical links in the combined field of the Sun and Earth. 

3.2. Flexing of the arms 

We also need to estimate the variation in armlength which is important for the TDI analysis to 
follow. Although there has been some adhoc work on this topic previously, the question needs 
to be revisited and a complete solution sought. Here we propose to address the question of 
the residual noise, having given the exact optical links. The first important task is to exactly 
compute the rate of change of arm-length. Figure |4] shows the rate of change of the six optical 
links as a function of time over a period of three years. 
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Figure 2. Comparison of flexing just due to Sun (dotted curve) and the combined effect of 
Sun and Earth (bold curve). The flexing of the arm is more and especially high in the third 
year in the combined field of the Sun and Earth as compared to that of Sun only. The figure 
is plotted for a duration of three years and with (f>o = for the link 12. 
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Figure 3. The figure shows the variation in the six optical links of the LISA model for three 
years. The lengths are given in metres. 

We find that in the optimised model of LISA configuration, this rate of change is less than 
4 m/sec. if we just consider the Sun's field. Including the Earth's field the flexing still remains 
< 6 m/sec in the first two years and increases to < 8 m/sec in the third year. Earlier estimates 
were ~ 10 m/sec. This numerical estimates are most crucial for their effect on residual laser 
frequency noise in the TDI. 

4. Time-delay interferometry with variable armlengths 

In order to cancel the laser frequency noise, time-delayed data streams are added together in 
which an appropriate set of time-delays are chosen. In general the time-delays are multiples 
of the photon transit time between pairs of spacecraft. In [4] a scheme based on modules over 
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Figure 4. The rate of change of armlengths for the six links is shown in units of m/sec. This 
rate of change is less than 6 m/sec upto the second year and increases to a maximum of about 
8 m/sec in the third year. 

commutative rings was given where the module of data combinations cancelling the laser noise was 
constructed. This fully cancels the laser frequency noise for stationary LISA. There are only three 
delay operators corresponding to the three armlengths and the time-delay operators commute. 
This scheme can be straight forwardly extended to moving LISA |8], where, now because of 
Sagnac effect, the up and down optical links have different armlengths (photon transit time) but 
the armlengths are still constant in time. Now there are six delay operators corresponding to 
the six optical links and they commute. These are the modified but still first generation TDI. It 
is crucial for the operators to commute if this scheme is to work - we must have a commutative 
ring. However, for LISA the armlengths do change as a function of time - flexing of the arms 
- and the first generation TDI modified or otherwise lead to imperfect cancellation of the laser 
frequency noise. Here we consider the model for LISA whose variation in arm-length is minimum 
when only the Sun's field is considered. This property may no more exactly hold in the combined 
field of the Sun and Earth, but assuming that Sun's effect is dominant the variation in armlength 
may still continue to be near optimal. Here, we however, assume the same initial conditions and 
compute the residual laser noise in the first generation modified TDI variables. 

Let C(t) = Av(t)/i/ represent the laser frequency noise in some optical link. Let Dj be the 
delay operator corresponding to the armlength Lj(t),i.e. DjC(t) = C(t — Lj(t)). If we operate 
on C(t) with operators Dj and Dk in different orders, it is easily seen that 

DjDkCit) + D k DjC{t). (21) 

The operators do not commute. A combinatorial approach has been adopted in [9] to deal 
with the totally noncommutative case. However, our aim here is to estimate the level of the 
noncommutativity of these operators in the context of the LISA model, compute the residual 
laser frequency noise and compare it with the other secondary noises in LISA for several of the 
TDI combinations. Our investigations will show whether the noncommutativity can be ignored 
or one must deal with non-commutative operators. Here our investigations necessarily take 
into account the symmetry of the LISA configuration and therefore we expect the residual laser 
frequency noise to be smaller than if the symmetries were absent. 
The first step is to develop a calculus of the Dj operators. 
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4-1. The calculus of the delay operators and relevant commutators 

We find that for this calculation we require to develop this only to the first order in L. This 
is because we find for our model L ~ 10~ 6 metres/sec 2 and thus even if one considers say 6 
successive optical paths, that is, about At ~ 100 seconds of light travel time, At 2 L <~ 10~ 2 
metres. This is well below few metres and thus can be neglected in the residual laser noise 
computation. Moreover, L 2 terms (and higher order) can be dropped since they are of the order 
of < 10 -15 (they come with a factor 1/c 2 ) which is much smaller than 1 part in 10 8 . The 
calculations which follow neglect these terms which simplifies the computations and makes them 
tractable. We begin with the effect of one operator on C{t) and then by induction obtain the 
effect n succesive operators operating on C(t). 

D k C(t) = C(t-L k (t))=E k C(t), (22) 

where E k is a delay operator with constant delay by the time L k (t) at the given time t. 

Applying the operators twice in succession and dropping higher order terms as explained 
above, 

Dk 2 D kl C = C(t — L kl (t — Lfe 2 ) — Lfe 2 ) , 

~ E k2 E kl C + L k2 L kl E k2 E kl C . (23) 
It is easy to generalise the above formula by induction to n operators: 
D kn ...D kl C = E kn ...E kl C + f n E kn ...E kl C , 

n p—1 
p=2 9=1 

It must be noted that both C and its time derivative are evaluated at the delayed time given by 
the successive application of n delay operators E km . 

We now turn to the commutators of the operators. These occur in many of the LISA 
observables and therefore it is useful compute these. The term in C cancels out; only the C term 
remains. The simplest of the commutators is: 

[Dj,D k ] = DjD k - D k Dj = L 3 L k - L k Lj , (25) 

where it is understood that the commutator multiplies C at the delayed time t — Lj (t) — L k (t) 
for fixed time t. 

For short we will write j instead of Dj whenever there is no possibility of confusion. Thus 
the commutator [Dj,D k ] will be simply written by \j,k]. We list few more commutators that 
occur in the observables: 

kmj - jkm = (L k + L m )Lj - L 3 (L k + L m ) , 

Imjk - jklm = (L t + L m )(Lj + L k ) - {Lj + L k )(U + L m ) , 

Imnxyz — xyzlmn = (Li + L m + L n ){L x + L y + L z ) (26) 

~ (L X + Ly + L z )(L l +L m + L n ). (27) 

This formula generalises in an obvious way to 2n operators. 

4-2. Residual laser frequency noise in some important TDI observables 

The laser frequency noise is usually assumed to be ~ 30Hz/\/Hz in most of the literature so far. 
However, by the time LISA flies the expectations are for this noise estimate to reduce to say 
Au <~ 10Hz/\/Hz. If we divide this number by the laser frequency f ~3x 10 14 Hz, we obtain 
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the noise estimate in the fractional Doppler shift C(t) = Ai/(t)/i/Q ~ 3 x 10 14 . Thus the power 
spectral density (PSD) of the noise C is: 

Sc(f) = (\C(f)\ 2 ) ~ KT 27 Hz" 1 , (28) 
where C(f) is the Fourier transform of C(t). 

From this equation it is easily deduced on differentiating that the PSD of the random variable 

C is: 

S d (f)=^ 2 fSc{f) Hz. (29) 
The commutators computed in the last subsection when divided by c 2 have dimension of time 
(assuming that the dot on L/. is just d/dt). A commutator is essentially the time difference in 
the transit times of photons along two different paths - the residual time. This is the reason 
why the laser frequency noise does not cancel out. If the two paths were exactly equal, the laser 
frequency noise would completely cancel out as it happens when the armlengths are constant. 
In order to compare the noise in a specific observable, say, the Sagnac or Michelson, one must 
compare the PSD of the secondary noise in that observable (which has dimensions of Hz -1 ) 
with the residual laser frequency noise PSD which is given by At 2 S^(f), where At is the time 
difference. Note that At is a function of t as the LISA constellation evolves and propagates in 
the gravitational field. We carry out these computations for the well known observables, such as 
the Sagnac, denoted by a,j3,j, the Michelson observable X, the symmetric Sagnac observable 
£. For the observable £, the commutators have different degree polynomials and therefore, the 
C appears at different delays. In the Fourier domain, these then appear as phase factors in At, 
which in effect becomes complex. 

We follow the notation and conventions of [8] and [4] which are the simplest for our purpose. 
The six links are denoted by U\V l ,i = 1,2,3. The time-delay for the link U 2 from S/C 1 to 
S/C 2 or 1 — ► 2 is denoted by x in [8] (which is 3' in [9] and so on in a cyclic fashion); the 
delay for link U 3 from 2 — ► 3 by y; the delay for link U 1 from 3 — ► 1 by z. The delays in 
the other sense are denoted by l,m,n. The delay for the link — V 1 from 2 — ► 1 by I; and then 
links V 2 ,V 3 and the corresponding delays m,n are defined through cyclic permutation. In the 
formalism any observable X is given by: 

X = p i V i + q i U i , (30) 
where pi, qt, i = 1, 2, 3 are polynomial vectors in the variables x, y, z, I, m, n. Thus X is specified 
by giving the six tuple polynomial vector (pi,qi). 

We observe the following approximate symmetries in our model: 

L x ~ Li, L y f=a L m , L z f=a L n , (31) 
which also implies (this combination occurs in the Sagnac observables), 

L x ~\~ Ly -\- L z ~ L\ L m -\- L n . (32) 

It was shown in [7] that only in the Sun's field cx (sin il(t — + k sin 30(t — where k 

is a constant and t\® are given constants. When we consider the sum L x + L y + L z , the phases 

for the links y and z namely, flt y °^ and flt z ^ differ by 27r/3 and 47r/3 from the phase of the link 
x and therefore their sum is close to zero. The same is true for the links /, m, n. Here, when 
we consider the combined field of the Sun and Earth, this is no longer true but we find that, 
\{L X + L y + L z ) — (Li + L rn + L n )\ < 1 m/sec and \L X — L;| < 0.8 m/sec upto the first three 
years in our model. The same is essentially true for the pairs of links y, m and z, n. Thus these 
pairs of operators essentially commute. 

Thus here, for practical purposes, we are not dealing with a set of totally noncommuting 
variables, but with an intermediate case in which the six variables partition pairwise into three 
pairs, such that for each pair the variables essentially commute, while any other combination of 
variables does not. 
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4-2.1. The Sagnac variables: The modified first generation TDI observable a is given by the 
polynomial vector in the form (pi, g,) by: 

a — (k, kI, Klm,rj,r]zy,r]z) , (33) 

where k = 1 — zyx and rj = 1 — Zmri. If the variables x,y, z,l,m,n commute then the laser 
frequency noise is fully cancelled. However, if they do not commute, there is a residual term. Let 
the laser frequency noises on each spacecraft i be C respectively (we consider a single effective 
laser frequency noise random variable on each spacecraft), then the residual term is: 

AC = a x Cx + a 2 C 2 + a 3 C 3 ■ (34) 

We find that a 2 — a 3 = and a% = [zyx, Imn] and so by Eq. |27|) : 

At (*) = ^2 t( L * + L v + L ^ L i + L m +L n )-(Li + L m + L n ) (L x + L y + L Z )}, (35) 

and thus AC = AtC\. Because the vary during the course of an year the At also varies 
during the year and so also the amplitude of the random variable AC. Thus the PSD of AC is: 

S AC (f; t) = 4TT 2 At(t) 2 fS c (f) ■ (36) 

This is the residual laser frequency noise in the observable a which depends on the epoch t. 

This noise must be compared with the secondary noise [2]. However, because we are 
considering the modified TDI Eq. (f33|) . there are extra factors k and r] which do not appear 
in the corresponding first generation TDI. These factors introduce an additional multiplicative 
factor, namely, 4sin 2 (37r/Lo) in the secondary noise PSD which leaves the SNR unchanged but 
must be considered when it is compared with the residual laser frequency noise given in Eq. ((361 . 
Thus, 

S a (f) = 4sin 2 (37r/L ){[8sin 2 3^/L + 16sin 2 7r/L ]5 acc + 65 opt } (37) 

where S acc = 2.5 x 10- 48 (//lHz)- 2 Hz~ 1 and S opt = 1.8 x 10~ 37 (//lHz) 2 Hz _1 . In the Figured 
we plot S a (f) and Sac(/;*) a t three epochs an year apart. 
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Figure 5. The 'top' curve shows the PSD S a (f) of the secondary noises. The straight lines 
are the PSDs of the residual noise at three epochs chosen an year apart. Clearly the residual 
laser noise is adequately below the secondary noises. 

We see that, clearly the residual laser frequency noise is an order or few orders of magnitude 
below the secondary noises depending on the epoch. Since the other Sagnac variables (3 and 7 
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Figure 6. The curve shows the PSD Sx{f) of the secondary noises. The straight lines are the 
PSDs of the residual noise at three epochs an year apart. The residual noises are essentially 
below the secondary noises. Also we note that in this case the straight lines lie very close to 
each other suggesting periodicity of a year. 

are obtained by cyclic permutations of the spacecraft, the residual laser noise is the same, if we 
assume that the noise PSDs of Cj are identical; the C\ is replaced by C 2 and C3 respectively 
for (3 and 7. Thus modified first generation TDI Sagnac variables suffice to suppress the laser 
frequency noise. 

As noted before, the basic reason for this remarkable cancellation is the symmetry inherent 
in the spacecraft configuration and physics as detailed in the equations Eq. (|3"Tj) and Eq. (|32|) . 

4-2.2. The Michelson variables: The TDI observable X is given by: 

X = (1 - zn, 0, (Ix - l)z, 1 - Ix, (zn - I);, 0) (38) 
The residual term is: 

AC = X X C X + X 2 C 2 + X 3 C 3 . (39) 
It is found that X 2 = X3 = and X x = [zn, Ix] and hence: 

At(t) = ^ + L n)( L i + U) - (Li + L X ){L Z + L n )\ (40) 

Using this value of At in Eq. (j36| and Eq. (f28|) gives the residual noise in the Michelson variable 
X. The PSD of the secondary noise is given by [2]: 

S x (f) = [&smHTTfL + 32sm 2 2TTfL }S acc + 16sm 2 27rfL S opt (41) 

The noise plots are shown in Figure [H Clearly at low frequencies / < 1 mHz the TDI variable 
X suffices, that is, the laser frequency noise is adequately suppressed. At higher frequencies in 
most of the frequency domain the residual noise remains below 20%. 

4-2.3. The symmetric Sagnac variables: The modified symmetric sagnac variables denoted by 
C in the literature split into three. We consider here just one of these because the noise in the 
others is essentially the same. It is given by: 

C = (y(zx — to), (In — y)z, (zx — m)l, m(ln — y), (In — y)z, (zx — to)/) . (42) 
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We denote it by just £, dropping the subscript, instead of £1, because we will not be explicitly 
discussing the other two cyclic permutations of £. The residual laser frequency noise term can 
again be written as AC = (1C1 + (2C2 + (3C3 where the £fc can be expressed in terms of the 
commutators, 

Ci = [m, y] + [In, zx] 
C, 2 = [{zx - m)l,y] 

C 3 = [m, {In - y)z] . (43) 

Note that here none of the £fc are zero and hence contribute to the total residual noise. Here 
it is more appropriate to compute the random variables ACk defined below in terms of the 
commutators given above and the delay operators E^: 

c 2 AC\ = [—L y L m + L m Ly\E m E y C\ 

+ \{L n + Li){L x + L z ) - [L x + L z ){L n + L{j\E x E z E n EiC\ 
c 2 AC2 = [L y {Li + L m ) — L y [Li + L m )]E m EiE y C2 

+ [L y (L x + L z + Li) — L y {L x + L z + Li)]E y EiE x E z C2 , 
c 2 ACs = [—L m (L y + L z ) + L m {L y + L z )]E m EyE z Cz 

+ [-L m (L n + L Z + L t ) + L m {L n + L Z + Li)]E m EiE n E z C 3 . (44) 

We note here that in each equation, the Ck is delayed by different amounts. In the Fourier space 
this is translated into complex phase factors. For the purpose of this computation assuming 
equal armlengths Lq for all the links, the effective complex Aii is given by: 

c 2 Ah{t) = [-L y L m + L m L v Y- AmSLo 

+ [{L n + Li){L x + L z ) - {L x + L z ){L n + U)] e -^ lfLa . (45) 

Similar expressions can be derived for Ai 2 and At 3 . When computing the PSDs it is the modulus 
of Atk that enters into their expressions accounting for the phases. 

We now assume that laser noises Ck are independent of each other and also have identical 
PSDs. Therefore we may add the noises quadratically - that is, we take the sum of the PSDs. 
The result is: 

S AC (f;t) = 4 7 r 2 / 2 (|Ai 1 | 2 + |At 2 | 2 + |Ai 3 | 2 )S c (/) . (46) 

This equation describes the residual laser noise. We must now compare this PSD with the sum 
of the PSDs of optical and acceleration noises. This PSD is given by: 

S c {f) = 4 sin 2 ^/L (24sin 2 irfL„S acc + 6S opt ) . (47) 

There is an extra factor of 4 sin 2 nfLo in the modified TDI variable zeta in the PSD which 
must be considered when comparing with the residual laser frequency noise. The noise PSDs are 
plotted in Figure [71 At the low frequency end / < 1 mHz, the residual noise is close to S^{f), 
while at higher frequencies it less by an order of magnitude and thus reasonably suppressed. 

5. Concluding remarks 

In this work we have included the effect of the Earth on the flexing of the arms of LISA. We have 
computed the spacecraft orbits in the combined field of the Sun and Earth approximately and 
from this deduced the flexing of the arms of LISA by choosing the model which gave minimum 
flexing when only the Sun's field was taken into account. We note that the flexing in the 
combined field is no more periodic as was the case when only the Sun's field was considered. We 
have ignored the effect of Jupiter because we believe this effect to be not so dominant as that 
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Figure 7. The bold curve shows the PSD S^(f) of the secondary noises. The rest of the 
curves are the PSDs of the residual noise at three epochs an year apart. The residual noises 
are below the secondary noises except near the low frequency end. 

of the Earth. Firstly, if we compute the tidal parameter ej for Jupiter, ej = GMj/dj, similar 
to e of Earth, where, Mj « 2 x 10 27 kg is the mass of Jupiter and dj, the distance from LISA 
to Jupiter, which we take on the average to be ~ 5 A. U., then ej/e ~ 0.09. Thus we expect 
the effect of Jupiter to be less than 10% than that due to Earth. Secondly, the forcing terms of 
Jupiter have the periodicity pertaining to its own orbit and therefore will not be in resonance as 
was the case with the Earth, so we do not expect the effect to accumulate in the first few years, 
when the perturbations are small. Note that these results are valid so long as we can neglect the 
nonlinearities arising from higher order terms in e and a. 

We have then used the results of the flexing of LISA's arms to compute the residual 
laser frequency noise in important TDI variables, namely, the Sagnac, the Michelson and the 
Symmetric Sagnac. Our results are obtained to the first order in L dropping terms of degree / order 
equal to or higher than L? and L. We have compared the residual noises with the corresponding 
secondary noises. We find that the residual laser noise in all these variables tends not to be 
very high. In the Sagnac variables it is negligible, in the Michelson variables it is less than 20%, 
while in C, variables only at the low frequency end the residual noise becomes comparable to 
the secondary noises. If this is acceptable then the modified first generation TDI observables 
could as well be used along with our model of LISA. Second generation TDI variables generally 
involving higher degree polynomials may not be then required. 

Our model of LISA is optimal (minimal flexing of arms) only in the Sun's field. Clearly this 
opens up the question of seeking an optimal model for the LISA configuration in the field of the 
Sun, Earth, Jupiter and other planets which will minimise the flexing of the arms and therefore 
the residual laser frequency noise in the modified first generation TDI. 

We finally remark that our computations here may be useful in the development of a LISA 
simulator. This is because our computation of the optical links have been carried out within a 
fully general relativistic framework and we have taken into account the gravitational field of the 
Earth. Also any discrepancy observed between actual data and the model may suggest physical 
causes which would be of interest and therefore would have to be incorporated into the data 
analysis. 
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